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Abstract We report an indication (3.22er) of « 1860 Hz quasi-periodic oscillations 
from a neutron star low-mass X-ray binary 4U 1636-536. If confirmed, this will be by 
far the highest frequency feature observed from an accreting neutron star system, and 
hence could be very useful to understand such systems. This plausible timing feature 
was observed simultaneously with lower (sa 585 Hz) and upper (« 904 Hz) kilohertz 
quasi-periodic oscillations. The two kilohertz quasi-periodic oscillation frequencies had 
the ratio of w 1.5, and the frequency of the alleged w 1860 Hz feature was close to the 
triple and the double of these frequencies. This can be useful to constrain the models of 
all the three features. In particular, the 1860 Hz feature could be (1) from a new and 
heretofore unknown class of quasi-periodic oscillations, or (2) the first observed overtone 
of lower or upper kilohertz quasi-periodic oscillations. Finally we note that, although the 
relatively low significance of the w 1860 Hz feature argues for caution, even a 3.22cr fea- 
ture at such a uniquely high frequency should be interesting enough to spur a systematic 
search in the archival data, as well as to scientifically motivate sufficiently large timing 
instruments for the next generation X-ray missions. 
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1 INTRODUCTION 

Quasi-periodic os cillations (QPOs) a re observed from many neutron star low-mass X-ray binary 
(LMXB) systems (Ivan der KlisL l2006h . There are several kinds of QPOs, such as millihertz QPO, hor- 
izontal branch QPO, normal/flarin g branch QPO, hectohertz QPO, lower kilohertz (kHz) QPO, upper 
kHz QPO, etc. (Ivan der Klisl 12006b . In this paper we will mostly concentrate on high frequency QPOs, 
that can currently be observed only with the proportional counter array (PCA) instrument of the Rossi 
X-ray Timing Explorer (RXTE) space mission. KHz QPOs are such high frequency QPOs, and often 
they appear as a pair. For a given source, their frequencies normally move up and down together in the 
f» 200 — 1200 Hz range in correlation with the source state. KHz QPOs are scientifically very important 
for the following reasons: (1) their high frequencies indicate that they originate from regions close to 
the neutron stars, and hence could be used as a tool to probe the strong gravity region around these 
stars, as well as the neutron star properties, and (2) they have been observed from many sources, and 
repeatedly from a given source, with high significance. Indeed, there are several models that involve 
relativistic or bital frequencies and the neutron star spin frequency, as well as beating and resonances 
among them dStella & VietriL Il998t iLamb et al. L fl985l iMiller et al.l 119981: lAbramowicz & Kluzniakl 
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2009). However, although there are many kHz QPO models available in the literature, none of them can 
explain all the major properties of this timing feature. Therefore, these QPOs cannot yet be used as a 
reliable tool to probe the strong gravity region, or to measure the neutron star parameters. 

An intriguing aspe ct of the kHz QPOs is that so far they have never been observed with a frequency 
greater than 1330 Hz (Ivan der Klisl 120061) . while PCA can easily detect a feat ure with a much highe r 
frequency. In fact, even this 1330 Hz QPO was not confirmed by a later analysis dBoutelier et al. ll2009l) . 
This implies that no timing feature has ever been observed with > 1330 Hz from an accreting neutron 
star system. In this paper, we r eport the indication of an w 1860 Hz QPO from the persistent neutron 
star LMXB 4U 1636-536 (see lAltamirano et al. I d2008l) for the other timing features). This plausible 
QPO was observed simultaneously with the pair of kHz QPOs, which suggests that the r* 1860 Hz 
QPO could be from a new and previously unknown class of very high freque n cy QPOs. Note that k Hz 
QPOs from 4U 1636-536 were first discovered by van der Klis et al. 1 dl996l): IZhang et al. I dl996allbb 
and later were reported by many authors ( Wiinands et al. L 1 99 7 |: | Va ughan e t al. L 11997c IZhang et al. , 
1997| :|Vaughan et al. lll998tlPsaltis et al.lll998l:lMendez et al. lll998clMendezLll998|lMisra & Shanthi . 
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2 DATA ANALYSIS AND RESULTS 



The LMXB 4U 1636-536 was observed with RXTE PCA on September 28, 2007. We analysed the data 
from the event file FS4f_19d86beO-19d8c749 (time: 16: 17:36 to 22:47:38, i.e., 23.402 ks duration) of the 
obsID 93091-01-01-000, and searched for QPOs. We excluded the data gaps and a burst, as well as the 
time intervals in which the number of used proportional counter units (PCU) was changed, and divided 
the 21.85 ks of data into M (= 2185) equal segments of 10 s duration. We performed a fast Fourier 
transform on each time segment for the entire PCA energy range, and the resulting power spectra (up 
to 2048 Hz) were averaged in order to reduce the noise. By combining W (= 4) consecutive frequency 
bins (making the frequency resolution 0.4 Hz), we found a tentative peak at the frequency 1861.45 
Hz (see the upper panel of Fig. 1). Before exploring this peak further, we binned the power spectrum 
more in order to search for kHz QPOs. This gave the clear evidence of two kHz QPOs at ks 585 Hz 
and « 904 Hz (see the lower panel of Fig. 1). For probing the properties of these QPOs, we fitted the 
power spectrum with a model, and minimized the corresponding \ 2 to get the best-fit parameter values. 
This model was the sum of a "constant" (to account for the Poisson noise level including the effect of 
deadtime), a "powerlaw" (to describe the red noise or broad-band feature at lower frequencies), and two 
Gaussians (to fit the kHz QPO pair). Then we divided the power spectrum with the continuum part of 
the best-fit model (see, for example, the lower panel of Fig. 1) and mult i plied by 2, in order to have 
the n oise distributed as x with 2MW degrees of freedom (Ivan der Klisl Il989l: iBhattacharvva et al. L 
I2006I) . The peak powers of the 585 Hz and 904 Hz QPOs were 2.00669 (for W = 1024) and 2.00959 
(for W = 2048) respectively, which gave the respective single trial significance of 2.87 x 10 -7 and 
2.20 x 10 - 24 . Therefore, multiplying by the number of trials (20 and 10), we got a significance of 
5.74 x 10~ 6 (i.e., 4.5er) and 2.20 x 10~ 23 respectively. These numbers show that both the kHz QPOs 
were very significant. The Gaussian modeling gave the RMS amplitudes ks 4.5% and « 7%, and the 
quality factors w 5.6 and w 7.3 for the lower and the upper kHz QPOs respectively. 

In order to estimate the significance of the observed peak at « 1860 Hz, we followed the same 
procedure as mentioned in the previous paragraph. After having the noise distributed as \ 2 with 2MW 
degrees of freedom, the peak power 2. 1094 gave the single trial significance of 2.47 x 10~ 7 . We consid- 
ered the number of trials to be equal to the number of frequency bins (= 5120) in the power spectrum. 
This is because, (1) we applied no cuts to the data (e.g., based on energy bands), apart from the standard 
filtering for burst and gaps; and (2) we did not search for QPOs in any other data set of 4U 1636-536. 
Note that we chose this particular data set in connection with a research on broad iron lines, and we 
happened to find the feature at r* 1860 Hz. Furthermore, since the frequency of high frequency QPOs 
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evolves with time (Ivan derKlisll2006l) . it is unlikely that a feature with the same frequency will appear 
in another observation. Therefore, it should be reasonable to consider only the current data set in order 
to estimate the number of trials. Hence, considering the number of trials = 5120, we found that the 
plausible QPO at « 1860 Hz was 3.22ct significant. The RMS amplitude of this feature was w 1%, and 
it had an extremely high quality factor (Q « 4650), as this QPO was observed with 0.4 Hz resolution. 
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Fig. 1 Power spectrum (averaged over 2185 spectra; 
each from a 10 s time segment) of the 21.85 ks RXTE 
PCA data from the source 4U 1636-536 (see text). 
Upper panel (power spectrum with 0.4 Hz frequency 
resolution): a peak near 1860 Hz is clearly seen. Lower 
panel (same power spectrum with 51.2 Hz frequency 
resolution): two kHz QPOs (near 585 Hz and 904 Hz) 
are clearly seen, ler error bars and a model (solid line; 
describing the continuum part of the power spectrum) 
are shown. The vertical dotted lines mark the frequency 
of the narrow w 1860 Hz peak of the upper panel, and 
its half and one-third frequencies. Note that the latter 
two pass through the upper and lower kHz QPOs re- 
spectively. 

In order to find out whether the plausible w 1860 Hz QPO was stronger in a part of the dataset, we 
examined the power spectra of nine unequal data segments. This division was naturally done in order 
to exclude the data gaps, burst, etc. (see the first paragraph of §[2}. We searched for a plausible timing 
feature at the known frequency of pb 1860 Hz in the power spectrum of each of these data segments. 
Although we found the indication of such a feature in several segments, it was strong in the power 
spectrum of only one segment (hereafter segment 1; time: 19:20:06 to 20:53:26, i.e., 5.6 ks duration). 
This power spectrum was calculated by averaging the spectra of M (= 560) equal segments of 10 s 
duration. By combining W (= 32) consecutive frequency bins (making the frequency resolution 3.2 
Hz), we found a peak at the frequency 1860.85 Hz (see the upper panel of Fig. |2). After having the 
noise distributed as x 2 with 2MW degrees of freedom (see the first paragraph of §[2), we got a peak 
power of 2.07. Therefore, for the search at a known frequency in this particular power spectrum, the 
« 1860 Hz QPO had the significance of 1.82 x 10~ 6 (i.e., 4.8ct). The RMS amplitude of this QPO was 
2.4%. Therefore, this QPO was stronger in segment 1 than in the whole data set. Here we note that even 
for the segment 1, the « 1860 Hz QPO occured simultaneously with the lower and the upper kHz QPOs 
(see the middle panel of Fig. |2). 
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Fig. 2 Power spectrum (averaged over 560 spectra; 
each from a 10 s time segment) of the 5.6 ks RXTE 
PCA data from the source 4U 1636-536 (see § EJ). 
Upper panel (power spectrum with 3.2 Hz frequency 
resolution): a peak near 1860 Hz is clearly seen. 
Middle panel (same power spectrum with 5 1 .2 Hz fre- 
quency resolution): two kHz QPOs (at the same fre- 
quencies shown in Fig. 1) are clearly seen. A tentative 
broad peak near 1860 Hz is also seen, lcr error bars 
and a model (solid line; describing the continuum part 
of the power spectrum) are shown. The vertical dot- 
ted lines mark the frequency of the narrow pa 1860 
Hz peak of the upper panel, and its half and one-third 
frequencies. Note that the latter two pass through the 
upper and lower kHz QPOs respectively. Lower panel 
(same power spectrum with 12.8 Hz frequency reso- 
lution): the narrow part (marked with a vertical dotted 
line) and the broad part (suggested from the powers 
of several adjacent frequency bins above the average 
power) of the plausible pa 1860 Hz QPO are visible. 



The very high Q-value of the plausible pa 1860 Hz QPO is somewhat unusual. We, therefore, 
searched for a broader feature at the known 1860 Hz frequency in the power spectrum of segment 1, and 
indeed we found a tentative peak with the centroid frequency of 1868.85 Hz (middle panel of Fig. |2). 
After having the noise distributed as \ 2 with 2MW degrees of freedom, the peak power 2.0109 gave a 
significance of 1.82 x 10 -3 (i.e., 3. lcr). This tentative broad feature had the RMS amplitude of pa 3.7%, 
and a more usual quality factor of sa 35. Note that, although such a broad feature is not significant in 
the power spectrum of the whole dataset, the lower panel of Fig. 1 does show a small broad peak at the 
frequency of the plausible narrow pa 1860 Hz QPO. Therefore, while the evidence of a broad feature is 
weak, we still explored the data a little more in order to see whether the plausible narrow QPO could 
be a part of a broad QPO. The upper and the middle panels of Fig. [2] exclusively show the plausible 
narrow QPO and the plausible broad QPO respectively. Therefore, in the lower panel of the same figure, 
we used an intermediate frequency resolution of 12.8 Hz, and found both the plausible narrow QPO 
(marked with a dotted vertical line) and the plausible broad QPO (suggested from the powers of several 
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adjacent bins above the average power). This indicates that the narrow QPO at « 1860 Hz could be a 
part of a broader QPO. 

Finally, the vertical dotted lines of the figures show that the frequency of the plausible « 1860 Hz 
QPO is close to the double of the upper kHz QPO frequency and triple of the lower kHz QPO frequency. 
Therefore, the w 1860 Hz QPO could be the first overtone of the upper kHz QPO or the second overtone 
of the lower kHz QPO. 



3 DISCUSSION 



In this paper, we report the detection of a pair of kHz QPOs from a neutron star LMXB 4U 1636-536. 
The ratio of these QPO frequencies is r oughly 1.5, whi ch is somewhat consis tent with t he resonance 
model described in Abramowicz et al.l d2003l) (but see iBoutelier et al. I d2009l) : see also IZhang et al. I 
d2006l) I It is generally observed that the separation of the twin kHz QPO frequencies clusters around 
the neutron star spin frequency, or half of that. Models have been proposed in order to explain this aspect. 
However, for the kHz QPOs reported in this paper, the separation is m 319 H z ; which is quite d ifferent 
from the half of the stellar spin frequency (see also lMendez & Bellonil (l2007h:iYin et al~l (120071)'). N ote 



that the neutron star of 4U 1636-536 spins with a rate of 582 Hz dStrohmaver & Markwardtll2002l) . 

Apart from the kHz QPOs, we have found an indication of a QPO at w 1860 Hz. This is by far the 
highest frequency feature observed from an accreting neutron star system, and hence could be extremely 
useful to understand such systems. Moreover, this plausible QPO was observed simultaneously with the 
lower and the upper kHz QPOs, which means that the w 1860 Hz QPO could be from a new class of 
QPOs. 

We will now briefly discuss the plausible origin of this QPO using general arguments, and 
without going into specifics. In some of the models, the lower and upper kHz QPO and some 
other QPO frequencies are thought t o be one or more o f the following frequencies (for equato- 
rial circular orbits in Kerr space time; Ivan der Klisl d2006l) ): (1) Keplerian orbital frequency v<j, = 
Vk(1 + jirg/rf/ 2 )' 1 = (2 ! r)- 1 (GM/r 3 ) 1/2 (l + i^g/r) 3 / 2 )" 1 ; (2) radial epicyclic frequency 
v r = — 6(r g /r) + 8j(r g /r) 3 / 2 — 3j 2 (r s /r) 2 ) 1 / 2 ; (3) vertical epicyclic frequency vg = i^(l — 
4j(r s /r) 3 / 2 + 3j 2 (r g /r) 2 ) 1 / 2 ; (4) periastron precession frequency ^ pcr i = — v r ; and (5) nodal 
precession frequency z^nodai = v^ — vg. Here the last four frequencies are for infinitesimally tilted and 
eccentric orbits, M is the neutron star mass, r g = GMj c 2 , j = Jc/GM 2 , J is the total angular momen- 
tum of the neutron star, and r is the distance of an orbit from the centre of the neutron star. For a sample 
of radio pulsars, M was found to be distributed around 1.35M in a narrow Gaussian (a = 0.04M Q ; 
iThorsett & Chakrabartvl dl999l) ). Therefore, considering Al = 1.35Af Q , we find that only v$ or vg can 
be the frequency of the plausible ps 1860 Hz QPO, and that too for relatively large angular momentum 
parameter values (as far as a neutron star is concerned) and for orbits very close to the innermost stable 
circular orbit (ISCO; see Fig. [3). As the neutron star mass increases due to accretion, M is expected 
to be greater than 1.35M© for accreting neutron stars (such as the one in 4U 1636-536). Moreover, 
since vk = (c 3 /2-irG)(r /r g )~ 3 / 2 (1/M), vk decreases as M increases for a given r/r g . Therefore, 
it is unlikely that any of the above mentioned five frequencies can be the frequency of the plausible 
« 1860 Hz QPO. Besides, since the spin frequency of the neutron star in 4U 1636-536 is f ap i n = 582 
Hz, the beating of 24 P i n with any of v$, v r or vg cannot explain the frequency of this QPO. However, 
the w 1860 Hz QPO could be an overtone of any of v$, v r , vg or f per i, or beating between an overtone 
and a fundamental of these frequencies. This plausible QPO could also be an overtone of one of the kHz 
QPOs (see §|2). 

Although it is unlikely that the r* 1860 Hz QPO frequency is a Keplerian frequency it is still 
instructive to examine what constraints = 1860 Hz can impose on the neutron star parameter val- 
ues. This is because, the measurement of neutron star parameters provides the only way to constrain 
the theoretically proposed equation of state ( EoS) models of neutron sta r cores, and hence to under- 
stand the nature of supranuclear core matter dLattimer & Prakashl 120071 : iBhattacharvva et a"Q l2000l 
120011) . Various authors suggested ways to const rain the neutron star parameters assuming one of the 
kHz QPO frequencies as a Keplerian frequency ( Kaaret et al. L 1 19971 : iMiller et al. L Il998t IZhang et al. L 
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120071: IZhanei l2009h . For example, this assumption, and th e following two rea sonable conditions can 
constrain the mass (M) and the radius (R) of a neutron star dMiller et al. Lll998l) : 

R < r, (1) 

where r is the radius of the orbit associated with the kHz QPO via the expression of v$\ and 

risco < r, (2) 

where nsco is the radius of the ISCO. This is because the first condition gives a mass-dependent 
upper limit on R via the expression of and the second condition gives an upper limit on M: Al < 
c 3 1 (2n6 3 / 2 Gi / c i > \ r ) (for Schwarzschild spacetime). Therefore, if 1860 Hz were a Keplerian frequency, 
then the upper limit of neutron star mass and radius would b e ~ 1.2M W and w 10.5 km respectively . 
Such upper limits would support the strange star EoS models dCheng et al. Ill 9981 : iBhattacharvva et al. L 
l200ll) . although the less exotic neutron star EoS models could not be completely ruled out. 




Fig. 3 Radial profiles of various frequencies (colour 
coded) of equatorial circular orbits in Kerr spacetime 
(see §[3]). Two angular momentum parameters (j = 0.0 
(solid) and j = 0.3 (dotted)), and neutron star mass 
= 1.35A/0 are used. Dashed horizontal line exhibits 
the 1860 Hz frequency. This figure shows that these 
frequencies, that are often used to explain kHz and 
some other QPOs, cannot possibly explain a w 1860 
Hz QPO. 



The uniquely high frequency of the plausible w 1860 Hz QPO and its coexistence with both the 
kHz QPOs make it extremely in t erestin g. Furthermore, although no such high frequency signal was 
claimed previously, Ivan der Ml (120061) mentioned, "the distinct impression of the observers is that 
there is still much hiding below the formal detection levels". It should, therefore, be worthwhile to 
report this feature, which will spur a systematic search of such high frequency features in the archival 
RXTE PCA data, as well as in the large area xenon proportional counters (LAXPC) instrument data 
of the upcoming Astrosat space mission. Moreover, next generation X-ray timing instruments, such as 
the high timing resolution spectrometer (HTRS; proposed for the International X-ray Observatory) and 
the proposed Si pixel detector of the Advanced X-ray Timing Array, will have much bet t er cap ability 
to detect weak QPOs (see, for example, Fig. 2 of iBarret et aU (120081) : IChakrabartv et al. I d2008l) ). The 
plausible w 1860 Hz QPO will therefore motivate the future X-ray timing instruments. 
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